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Quasi-periodic Oscillations in XTE J1550-564: The rms-flux relation 
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ABSTRACT 

We present an analysis of the short timescale variations in the properties of the strong (type 
"C") quasi-periodic oscillation observed in XTE J 1550-564 during its 1998 outburst. In par- 
ticular, the QPO shows a correlation between absolute rms amplitude and mean source flux 
over timescales shorter than ~ 3 ksec. A linear rms-flux relation has been observed to be a 
common property of broad-band noise but here we report the first detection of rms-flux de- 
pendence in a QPO. The gradient of the rms-flux relation is correlated with the QPO peak 
frequency: from a strong positive correlation when the QPO peak frequency is below ~ 4 
Hz, through no correlation, to a strong negative correlation when the peak frequency is above 
6 Hz. This is the first time a negative short term rms-flux relation has been observed in any 
component of the power spectrum. Previous work on both the broad-band noise and QPOs 
in a range of sources have suggested the presence of a filter reducing the amplitude of QPOs 
with increasing frequency. We attempt to remove the possible effects of this filter and find that 
the previously negative rms-flux relations above ~ 5 Hz become constant. 
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1 INTRODUCTION 

Accreting black holes show variations over a broad range of 
timescales. In the well-studied Galactic Black hole binaries 
(BHBs) their power spectra may include broad-band noise, band- 
limited noise (BLN) and quasi-periodic oscillations (QPOs), dis- 
tingu is hed by their relative frequency w idths (see e.g. Ivan der Klis 
l2006t IMcClintock & Remillard l200d for reviews). The exact 
origin of the variability and its characteristic frequencies re- 
main elusive, although strong correlations between frequencies 
of different power-spectral components have been observed as 
well as corr elations between frequencies and spectral parame- 
ters (see e.g.lKlein-Wolt & van der Klisll2008l : IBelloni et al.1l2002l : 
iRemillard et alj|2002h . 

QPOs are observed over a range of frequencies in BHBs and 
NS, in BHBs they are commonly desc ribed as either low frequency 
< 50 Hz or high frequency > 100 Hz ( Ivan der KlisfeoOol) . Low fre- 
quency QPOs are often observed with further harmonics and sub 
harmonics and their formation mechanism is un clear. The three 
main types are commonly classified as A, B o r C (Wijnands et al. 
19991 : iHoman et al.ll200l[ ; IBelloni et alj|2002l ; ICasella et alj|2004 
20051) . with types A and B observed in the soft intermediate states, 
and displaying different harmonic and phase lag behaviours, and 
the type Cs observed in hard states (i.e low-hard and hard interme- 
diate) with strong broad-band noise. The frequency of type C QPOs 
is known to correlate with spectral properties of the source, the 
photon index, disc inner radiu s, disc temperature and both power 
law a nd disc fluxes (see e.g. [Remillard et al. 2002; Mun o et ail 
1 19991) although there is as yet no consensus on the physical ori- 



gin of QPOs dvan der Klisll2006l : IMcClintock & Remillard 1 20061 : 
iDone et alj2007l) . 

One of the simplest and most stable relationships observed 
between variability properties is the rms-flux relation, which con- 
nects the rms amplitude of variations to the mean flux level 
by a positive linear correlation that appears to operate over a 
wide range of timescales where the power spectra remain sim- 
ilar JUttlev & McHardvl 1200 ll : Ideissner et ail |2004) . This rela- 
tion is inherent to the short-term variability of sources with an 
otherwise-stationary (or close to stationary) power-spectral shape. 
Longer-term changes in rms can be caused by changes in power- 
spectral shape which correlate with the spectral evolution of the 
source within and between the different spectral states, but these 
are different to the rms-flux relation which appears to be a fun- 
damental aspect of the variability process itself. The linear rms- 
flux relation was initially observed in the BHB Cygnus X-l, 
to date the same linear rms-flux relation has been observed in 
Active Galactic Nuc l ei duttlev & McHardvl l200ll ; IVaughan etail 
120031 : lGaskellll2004l: lUttlev et alj|2005h. one n eutron star X-ray 
binary JUttlev & McHardvl 1200 ll; luttlevl l2004h and one Ultra- 
Luminous X-ray source (ULX) JHeil & Vau ghan 2010), it has also 
recently been observed ove r long timescales in the BHB GX 339-4 
( Mun oz43arias et al] [20101. This relation is observed in the noise 
components of the power spectrum, and in t he 401 Hz pu lsation 
of the neutron star system SAX 11808.4-3658 JUttley||2004l) . but to 
date has not been studied for QPOs. 

Arguably the most promising explanation for the linear rms- 
flux relation observed for the broad-band noise of accreting sys- 
tems is the 'propagating fluctuation' model (see discussion and 
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references in lUttlev & McHardvl f200lh : iGleissner et"afl d2004l) : 
lUttlevI d2004h ; lUttlev et al.l d2005h ). In this general model long 
timescale fluctuations in the accretion rate originate from large 
radii in an accretion disc (perhaps as random variations in viscos- 
ity) and propagate inwards, mo dulating shorter timescale variations 
origina ting at smaller ra dii (e.g. lLvubarskiilll997l:lKing et alj2004 
lArevalo & Uttlevll2006h . and eventually modulating the accretion 
rate in the innermost, X-ray producing regions. This simple scheme 
naturally explains X-ray emission that varies on a wide range of 
timescales, and a linear rms-flux relation due to the multiplicative 
coupling of variations on all timescales. It is not clear how QPOs 
fit into this scheme. 

Secular (~ 1 day) variations in the QPO properties (rms and 
frequency) are known to be correlated with the (energy spectral) 
evolution of black hole transient outbursts (i.e. changes in the 
physical state of the system). But this overall evolution of the 
outburst properties tends to follow distinct patterns, such as the 
hysteresis curve on the hardness- intensity diagram (HID; see e.g. 
iBelloni et al] 120051 : iBellonil |2010|) . As previously discussed, the 
within-observation variations (timescales < 3ks studied here) are 
variations in an almost stationary process (presumably random fluc- 
tuations within the accretion flow), the average properti es of which 
evolve through outburst (see lMunoz-Darias et alj|201ol . and Figure 
[T}. It is these relatively small and flux-correlated deviations from 
stationarity, and how they relate to other source properties, that are 
the subject of this paper. 

We focus on the type "C" QPOs observed at the begin- 
ning of the the 1998 outburst of XTE J1550-564. The long- 
ter m properties of this source ha v e bee n studied extensivel y 
bvlCui et alj dl999h.lWiinands et all f l999). Ho man et alj |200l[) . 



Remillard et al. (2002), Sobczak et al. (2000b), Reillv et al. (2001), 
Chakrabarti et alj d20091 ). and IRao et alj d20ld) . ISobczak et alj 
(2000b) compared the changes in amplitude and frequency for the 



QPOs over the course of the outburst to the energy spectral com- 
ponents in both XTE J1550-564 and the 1996-97 outburst of GRO 
J 1655-40. They found that the QPO frequency appears to be linked 
to both the disc flux and the strength of the power law. They also ob- 
serve that for the outburst of XTE J 1550-564, the amplitude of the 
QPO is closely linked to the QPO frequency: it rises until around 
3 Hz but appears to then fall off above this point. IRao et all d2010l) 
split each observation into 128s segments and demonstrated that 
over time the frequency of the fundamental QPO component can 
be linked to the disc count rate suggesting that there is some cor- 
relation between frequency, flux and QPO amplitude over short 
timescales. We further examine the short term variations in these 
parameters for the QPO, concentrating on the rms-flux relation 
within each observation and its relation to the peak frequency. 



2 OBSERVATIONS 

We use observations taken by the PCA detector on RXTE from the 
1998 outburst of XTE J1550-564 in programs P30188 and P30191, 
covering the first 46 days of the outburst. The data modes used all 
had 8 bit wide counters and are not likely to suffer the buffer over- 
flow problems which can dis tort the rms-flux relat ion away from 
a linear relationship (seen in IGleissner et alj|2004t Appendix A). 
During the first half of the outburst the source reached a very bright 
intermediate state where the strong ("type C") QPO shows high 
rms amplitudes, this means that they are ideal candidates to test the 
behaviour of the rms-flux relation for QPOs and how this relates to 
the continuum. 



Binned mode data has been used for this analysis. There are 
a few variations of this data mode used over the outburst, so the 
energy ranges for some of the observations differ. The ranges have 
been chosen so that they are as consistent across the sample as pos- 
sible, although slight variations remain between observations. The 
energy range chosen was as close to 4-14 keV as possible, and the 
difference between energy bands is no greater than 0.5 keV. 



3 DATA ANALYSIS 

In order to measure the rms-flux relation each continuous light 
curve is divided into non-overlapping 3 s intervals. Due to the lower 
frequency limit imposed by this segment size we only use observa- 
tions from this sample where the primary type C QPO is above 1 
Hz. In practice this includes most of the observations where the 
QPO is detected, only excluding six, leaving 48 observations to be 
processed. From each 3 s segment the power spectrum was esti- 
mated in absolute normalisation using standard techniques, and the 
background subtracted mean count rate of the interval, corrected to 
1 PCU detector, was recorded. The power spectral estimates were 
then averaged in non-overlapping flux bins, with at least 20 inter- 
vals contributing to each bin, to provide well-determined power 
spectra as a function of flu x. These power spectra were then fitted 

in xspec 12 dArnaudl 19961) . 

The models discussed in IRao et alj JioToh have been applied, 
whe re multiple Lore ntzians are used to describe the fit (see figure 
1 of IRao et alj|2010l for representative fitted power spectra). In all 
observations the model used consisted of four Lorentzians: Two of 
these were fitted to the broad band noise, they are referred to in 
IRao et all |2010) as L ft and L„ n (ft stands for flat topped noise at 
low frequencies < 1 Hz and pn for peaked noise, generally seen 
at frequencies close to, or just above, the QPO harmonic). These 
two Lorentzians had Q-factors less than 2{Q = v/Av where v is the 
cen troid frequency and Av the full width at half maximum). Follow- 
ing lRao et alj d2010h the primary QPO component is defined as the 
strongest feature in each observation. Due to the limited frequency 
range (0.33 - 100 Hz) within the flux-binned power spectra it is 
rarely possible to resolve the fundamental and both the harmonic 
and subharmonic simultaneously. For this reason the two further 
Lorentzians described the primary QPO component and either the 
harmonic or the sub-harmonic, depending on which was most visi- 
bl e. These Lorentz ians are referred to as L F and either L s or L h . As 
in lRao et al. (2010) the Poisson noise component was not removed 
from the power spectra, but is included as an additional constant in 
the model. This model provided an adequate fit to all observations 
within the sample, the reduced chi-squared varied from 0.75 to 
1.41, with the wo rst fit having a ch i-squared of 354.1 with 179 de- 
grees of freedom. IRao et all d2010l) indicated the need for a further 
high frequency broad-band component in some observations within 
their sample. All of these observations are ones where the QPO is 
below 1 Hz, therefore this was not necessary for the fitting of any 
observations used here. For each fit the Lorentzian properties were 
recorded, namely the total rms (in absolu te units, i.e ct/s/PCU), the 
peak frequency (v F ) and the width. See IPottschmidt et alj d2003h 
section 3 for appropriate formulae. The rms value is corrected to 
represent the power in the QPO over only positive frequencies. This 
normalisation is thus the integrated power in the QPO (R 2 F ) and the 

measured rms within the component is therefore a = yj^- 
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Figure 1. a. Rms-flux relations for the fundamental QPOs in each obser- 
vation within the sample. Points within observations are joined by dotted 
lines. Some observations with rms values close to 350 ct/s/PCU and mean 
count rates between 2000-2500 ct/s have been removed for clarity, b. Mea- 
sured fundamental QPO frequency for each of the flux bins used for the 
QPO rms-flux relations. The steeper change in frequency with flux at high 
frequencies is clearly visible. For both a. and b. the colour bar describes the 
average frequency of the fundamental QPO component. 



4 RESULTS 

Fig. Q] shows the dependence of the QPO amplitude (rms) and fre- 
quency on flux (count rate) both within and between observations 
(typically one day apart). Clearly the QPO frequency is varying 
on short timescales and is well corre l ated with the flux as m en- 
tioned previously bylRao et alj feoid) . ISobczak et ail fcOOObl) and 
lRemillardetZlj2002l) . the gradient of this relation appears to in- 
crease with average frequency. The rms shows a more complex re- 
lation with the flux: when the QPO is below 4 Hz it follows the pos- 
itive linear rms-flux relation observed in the broad-band noise, be- 
tween 4.0-5.5 Hz the rms becomes constant with flux, above 5.5 Hz 
the rms-flux relation is negative. This is the first time a short-term 
negative rms-flux relation has been observed. The average QPO fre- 
quency for an observation is not monotonically related to either the 
mean rms or flux over long timescales, even though clear relation- 
ships between these variables are found within each observation 
(see Figure 1). This illustrates the different effects of the long term 
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Figure 2. Gradient (k) of the rms-flux relation against frequency of the fun- 
damental QPO measured in both the fundamental (diamonds - black) and 
the harmonic (crosses - red). Both clearly follow a very similar pattern. 



evolution of the system vs. short term random variations in the ac- 
cretion flow. 

Fitting the rms-flux relations for each observation with a lin- 
ear function gives a reasonable measure of the gradient (k F ) which 
is clearly anti-correlated with the average QPO frequency (v F as 
seen in Figure [2] The point where the gradient reverses from posi- 
tive to negative occurs around 5.5 Hz. The gradient of the rms-flux 
relation for the harmonic, k h , changes simultaneously with that of 
the fundamental, k F . Figure 2 shows k/, becomes negative when k F 
does, i.e. around v F > 5.5 Hz (where v h > 11 Hz) (Figure [2] red 
points). We have been unable to asses the flux dependence of the 
high frequency power spectral continuum due to uncertainties in 
subtracting the power of the QPO and harmonics. 



5 DISCUSSION 
5.1 Summary 

We have shown that the rms amplitude (o>) and peak frequency 
(v F ) of the "type C" QPO of the microquasar XTE J1550-564 
(observed during the 1998 outburst) both vary strongly on short 
timescales (3 - 3000 s) and are correlated with the source flux. The 
rms-flux relation for the QPO depends on the peak frequency of 
the QPO: when v F ~ 1 - 5 Hz there is a positive, linear correlation 
between the QPO rms and flux, at around v F ~ 5 Hz the rms-flux 
relation becomes flat (rms independent of flux), and at higher fre- 
quencies the rms and flux becomes negatively correlated. For all 
values of v F examined (~ 1.5 - 10 Hz) the rms-flux relation is 
approximately linear, and the gradient itself appears to be mono- 
tonically related to the QPO peak frequency. Above ~ 3 Hz v F also 
increases with flux in an approximately linear manner, below this 
point it remains constant. The correlations between the flux and the 
QPO rms and peak frequency on these short timescales are distinct 
from the long term evolution of the QPO during the outburst that is 
apparent between observations separated typically by ~ 1 day. 

Within the framework of the "propagating fluctuation" model 
discussed in the Introduction, it is reasonable to imagine the QPOs 
as an enhancement in the variability amplitude at certain frequen- 
cies (radii) but still coupled to the broad-band spectrum of varia- 
tions. This would naturally lead to a linear, positive correlation be- 
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Figure 3. a ) Rms-flux relations for the QPO following the removal of a filter 
reducing the amplitude of the QPO with increasing frequency, b) Gradients 
found from fitting the above rms-flux relation with a linear function. Fre- 
quency dependence of the rms-flux relation is still visible, but the gradients 
no longer become negative. 



tween long term flux (driven by variations at frequencies lower than 
the QPO) and QPO rms amplitude. It is less obvious what causes 
the frequency-dependent changes in the rms-flux behaviour of the 
QPO. We explore some possibilities and implications below. 



5.2 The effect of a Frequency Dependent Filter 

Using these observations ISobczak et al. (2000b) showed that on 
long timescales the fractional rms of the QPO rises with frequency 
until around 3 Hz, above this it f alls off steeply (see t heir fig- 
ure 6). |Pottschmidt et aljf2003h and lAxelsson et al.|j20o3) showed 
that the broad Lorentzian components of the power spectrum of 
Cygnus X- l also show a sharp drop in amplitude above 3-5 Hz (see 
figure 6 of IPottschmidt et al]|2003l) . Similar correlations between 
QPO amplitude and frequency have been observed in the BHBs 
H1743-322, GRO J1650-500, XTE J1655-40 and GRS 1915+105 



jMcClintock et al]l2009l : iKalemci et ai1l2006l. 120031 : iDebnath et ail 
l2009l ; lMunoetaljl999l) . 

One explanation for the attenuation of these components as 
they move to higher frequencies is in terms o f a low pass filter 
acting to suppress va riability above ~ few Hz dPone et alj|2007l ; 
iGierlinski et al.l r2008). Proposed physical mechanisms for the sup- 
pression of fast variations include the damping of accretion rate 
perturbations faster than the viscous timescale in the inner disc 



dPsaltis & Normanll2000l ; lDone et alj|2007t) and scattering of emit- 
ted X-rays by an optic ally thick outflow created at high luminosity 
iTitarchuk et ail (2007) . Such a filter would affect the observed rms- 
flux behaviour for QPOs that move in frequency above ~ 3 - 5 Hz. 
Is the change in k F simply the result of a QPO frequency that in- 
creases with flux on short timescales, combined with the rms being 
more strongly suppressed (by the action of the filter) at higher fre- 
quencies and thus fluxes? In order to investigate this possibility we 
modelled the decay of the QPO (fractional) rms with frequency in 
terms of a doubly broken power law, and use this to "recover" the 
unfiltered QPO rms. 

A doubly broken power law over 1.5-10 Hz, with indices of a 
= -0.1, -0.5 and -1.7 and breaking at 3.4 Hz and 5.1 Hz respectively, 
can be used to describe the trend in the long term rms-frequency 
behaviour of the QPO. We note that this is likely an over-estimate of 
the effect of a filter, since it predicts filtering of the power spectrum 
by a power law with with index -3.4 above 5 Hz; the observed 
power spectrum is usually flatter. Assuming this function describes 
the filter it is possible to "correct" each data point in Figure Q] for 
the rms-reducing effect of the frequency-dependent filter. 

The rms-flux relations using the "filter-corrected" amplitudes 
are shown in Figure [3] There are two significant differences be- 
tween this and the raw data shown in Figure QJ: at higher fre- 
quencies the rms-flux gradient tends to zero above ~ 5 Hz, rather 
than becoming negative, and the short timescale rms-flux data 
(within each observation) now follow a similar curve to the longer 
timescale (between observation) changes in rms and flux. Signifi- 
cantly for the present discussion, the rms-flux relation still shows a 
monotonic dependence of gradient on the QPO frequency. A flat- 
tening to zero gradient might indicate a ceiling on the QPO strength 
(in absolute units) which effectively prevents the rms increasing 
past ~ 600 ct s when v F > 5 Hz. 

5.3 Physical implications of QPO behaviour 

Although we still lac k a complete ph ysical picture for the ori- 
gin of QPOs (see e.g. Ivan der Klisll2006l) . recently proposed ideas 
incl ude Lense-Thirring precession of an extended, hot flow (see 
e.g. [Ingram et al.l|2009h and a magnetohydrodynamic dynamo cy- 
cle dO'Neill et alJuOlQ i). However, it is not clear what, if any, pre- 
dictions such models make for the rms-flux relation of the QPO. 
One would hope that as these and similar models are developed 
and explored they will yield specific predictions about the coupling 
of QPO amplitude and the lower frequency noise variations that can 
then be compared against the observed behaviour. 

We have shown there is a clear change in the flux-dependence 
of the QPO that is itself linked to the QPO frequency, and that 
peak frequencies of v F « 5 - 6 Hz appear to stand out as mark- 
ing where the rms-flux gradient becomes zero. (This seems quite 
robust to the possible effects of a frequency-dependent filter.) This 
frequency range is already noteworthy as the location of the type B 
QPOs in this outburst and others. Given this coincidence it seems 
reasonable to speculate that 5-6 Hz may represent a relatively sta- 
ble characteristic frequency of the system (in contrast to the vari- 
able peak frequency of type C QPOs). Replication of these results 
for other sources and outbursts is necessary to establish whether a 
fixed characteristic frequency can be identified in this way. 



5.4 Correlations with the energy spectra 

It is well known that the properties of the type C QPOs in BHB 
outbursts (notably v F ) are strongly correlated with the various en- 
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ergy spectral param eters over lone times c ales, i.e. between ob- 
servations (see e.g. ISobczak et al. 2000b ; iRemillardetaLl |2002| ; 
iRao et alfcOld. for discussion of the 1998 outburst ofXTE J1550- 
564). U sing the spectral param eters obtained by both lSobczak et all 
( 2000a) and lDunn et all hold) we find strong correlations between 
the energy spectrum (e.g. power law index, disc flux, etc.) and the 
gradient k F of the rms-flux relation for the QPO, as would be ex- 
pected given that all these variables are strongly correlated with v F . 
At the present time the physical basis of these relations remains 

unclear. 

iDunn et alj d20ld) noted that this outburst of XTE J1550-564 
is unusual in that it reaches a particularly high intermediate state, 
at the peak of the outburst they find that the luminosity is close 
to l^Edd- It is during these high flux observations, where the QPO 
reaches its highest frequency, that the loss in positive correlation 
between rms and flux is most clearly seen. Further investigations 
into the effect of high luminosities, and their associated states, on 
the short-timescale behaviour of low frequency QPOs in other out- 
bursts and sources would be necessary to fully understand any link. 
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